Two approaches for Helmholtz equation: generalized 
Darboux Transformation and the method of 9-problem 



E.Sh.Gutshabash 

o . 

O ■ Institute of Physics, St. -Petersburg State University, Russia 

■ e-mail: gutshab@EG2097.spb.edu 

a ■ 

■ Abstract 

, Two approaches to solution of the two-dimensional Helmholtz equation with a "wave 

CO ' number" are proposed. The results can be applied both in numerical areas of physics and 

in the theory of nonlinear equations. The first approach is based on the requirement of 
the covariance of equation under the generalized Darboux transformation (Moutard trans- 
formation) . It allows to construct a new solution of equation, using a given initial solution 
' of the equation. Simultaneously we obtain the "dressing" relation for the "wave number". 

The simplest examples of the approach are considered in detail. In the second approach the 
Green-Cauchy formula (the 9-method) is applied to reduce the solution of the equation to 
the solution of a system of singular integral equations. 

> 

^ ■ 1. INTRODUCTION 

■ An intensive development of the methods in the theory of nonlinear integrable equ; 
, stimulated as well the development of the so-called dressing methods for the solution of a 

I of linear equations having the important physical applications. One of the most popula 

effective methods is the method of Darboux Transformation (DT) which gives the 
solutions of different one-dimensional equations. The situation for the multidimensional 
pH I tions is more complicate. In [2j the multidimensional variant of DT, which depends on s 

• derivatives, was considered, but such approach is rather complicate and cumbersome. 

In paper |2| the method of integration of the non-stationary Schrodinger and Fokker-^ 



X 



^ , equations was proposed with using of the Darboux-like anzats and introducing of some funct 

' differential equation. It appeared, that in some sense this approach is similar to the so-- 

Moutard transformation (see, for example, ||4j) for the variable real "wave number" (potei 
In this case the functional-differential equation of the second order is replaced by the syst 
two differential equations of the first order for functional, included in Darboux's anzats. F( 
complex "wave number" the situation is more complicate - the solution can be found in 
simple cases only. 

In Section 2 of this paper the idea of a Darboux-type anzats - Moutard transformatic 
introduced in the example of the two-dimensional Helmholtz equation. The main advant; 
such approach is independence of the procedure on the spatial derivatives, and hence it c 
straightforwardly extended onto higher dimensionality of space. 

In Section 3 the simplest classical problem of diffraction of the electromagnetic wa\ 
a half-plane is considered. We reformulate this problem in terms of a so-called meth 
9-problem. The idea of such approach (in another formulation of the problem) belor 
V.D.Lipovskii 0. The powerful and effective possibilities of this method for integration oi 
tidimensional systems such as Kadomtsev-Petviashvili-2 equation, Devey-Stewartson-2 s; 
and so on were demonstrated (see, for example, [Hl-IZI)- Apparantely, up to now this app 
was not applied in the theory of diffraction. Similar to the Wiener-Hopf method, it alio 
formulate the initial boundary problem in terms of the integral equations, which admit tl 



2. The generalized Darboux Transformation 

a) . The two-dimensional Helmholtz's equation belongs to a wide class of the equatic 
the following form (see, for example, jH]): 

where V = V{p, q) is the given, in general case, complex-valued function. Here and bek 
use the notations: = d"^/dq, .... 

In the simplest particular case: p = x, q = y, V{x,y) = 1, this equation has the g( 
solution 

'^ix,y) = / fiis)Joi2i^yix-s))ds+ /' f2is)Joi2iy/x{y - s) )ds+ 
Jo Jo 

+[/i(o) + /2(o)]Jo(2iViy), 

where Jo(-) is the Bessel's function, and /i, /2 are some differentiable functions. 

The class of such equations is important both for the differential geometry and for a 
variety of the physical problems. Indeed, for example, for p = x + t, q = x — t we ha\ 
heterogeneous wave equation: 

^tt - "^xx = V{x,t)'if, 
and for p = z = x + iy, q = p we obtain the Helmholtz equation itself: 

A^' = V{x,y)'^. 

We shall use it also in the complex form {dz = {l/2){dx — idy), dz = {\/2){dx + idy)) 

^z-z = \v{z,z)-i'. 

The equation (0^5)) arises in the Quantum Mechanics (two-dimensional stationary Schrc 
equation), theory of diffraction of electromagnetic waves, an acoustic diffraction and other 
lems. From the nonlinear equations point of view the interest to @ is caused by his n 
the associated linear system for the well-known nonlinear completely integrable Novikov-V( 
equation [Sj. 

It is not difficult to check by direct calculation that if $ = ^{z, z), x = xi^, z) ar 
solutions of equation (@J), then two important relations are true: 

{^zX)-z = {^-zX)z, {^zX - ^Xz)z = {^Xz - ^zX)z- 

It follows from the second of them, that the integral 

u;($, x) = / (^^X - Xz^)dz + {xz^ - X^z)dz, 

does not depend on the integration path. 

b) . We shall proceed to the analysis of the equation Q from the generalized Dai 
Transformation point of view. Let are two linearly independent solutions of ©. Folli 
^ |3j , we assume that 



where = ^i) = 0,(z, z) is the functional, given on direct product of two copies of s 
of "wave functions ". After substitution of the anzats in equation 0, and requirin 
covariance of this equation under transformation ^' —> [/—>[/, we obtain the equ 
which do not contain already the potential: 

- {in^i)sn, - {in^i),n2 + gi^ = o. 

Here qi = qi{z, z) is a complex- valued function ( "parameter" of separation of variables), 
duced in order to obtain the real potential V{z, z) for the case of V{z, z) = V{z, z) (conser^ 
medium). If V{z,z) ^ V{z,z) and V{z,z) ^ V{z,z) (nonconservative medium), one maj 
gi = in 0. 

Equation @ allows the physical interpretation: it looks like the two-dimensional stati' 
Schrodinger equation for charged particle in a non-homogenious stationary electromagneti( 
^. Actually, this equation can be written as (see, for example, [TO]): 

AV- - 2i(AVV') - (A^ + - E)ip = 0, 

where A(x, y) = {Ai, A2) is the vector potential, with the gauge condition divA = 0, cj) = (j) 
is the scalar potential, E is the energy value, and the system of units h = c = e = 1, m - 
was used (m is the mass of the particle). In variables z, z : 

i^z-z - \{iAi - A2)^, - ^{iAi + A2)^Ps - ^[A^ + (I) - E]^ = 0. 

Comparing (fTT|) and @, we see, that analogues of coefficients at Qz, ^2 are expressed as 
combinations of the components of the vector-potential, and the analogue of function q 
terms of the square of its module and of the scalar potential. 

Besides, the separation of variables leads to the following dressing relation for potential 
number) : 

Viz,z)=Viz,z)-8{ln^fi)z,-iqi{z,z). 
In the conservative case the requirement of reality reads: 



{arg '^i)zz = -4:1m qi. 
For the double dressing procedure, with 

and taking into account (|12j) . we obtain 

V{z, z) = V -8{lnQ{^u ^2)).s-4{gi + g2[l]}, 

where ^2[1] is the fixed solution of equation @ with ^ ^[1]; ^ ~^ ^[1]) and ^'2 
^2)/^i7 ^2 is some fixed solution of equation Q. 
It is clear, that similar to the case of standard DT, the dressing procedure proposec 
can be multiply iterated. Setting = ^([0], ^ = -^[l], . . . ,V = V[0], = . . ., we ha^ 
chains: ^-[0] ^[1] ^ . . . ^ *[7V] and V[0] V[l] . . . ^ V[N] .... In parti 

for the A^-fold iteration of the dressing procedure, one will obtain the relations, exprese 
terms of the functions ^(^i, ^i) only [T]. 



c) . Now some simplest examples of realization of the proposed approach will be consic 
1. Let us take V{z, z) = in ©. Then ^j{z, z) = 'iff\z) + ^f\z), \ '^f^ are ab 
functions, j = 1, 2, . . . , A^, . . .. The relation (|12() takes a form: 

V[l] = -8[ln (^/f + ^?{z))U - 
In the case V[l] = V[l] the following condition 



Im qi 



(^(l)+^(2))2 ($(!)+ $(2))2 

must be satisfied. On the next step we have: 

V[2] =V- 8{lnn{^i, ^2))zz - 4[gi + q2[l]), 

and the functional 'I'2) satisfies the equation ©, if to replace ^' ^'i, — > ^^21 Qi ^ 

in it. It is not difficult to check, that for qi = q2[^] = this equation is solved by the fui 
^'2) = ^2)) where the expression for <^(^'i, ^'2) is defined in After the s 

calculations we obtain: 

n{^i, ^2) = 70 + 2(^;^)M/f ) - ^'^^^f V 

where 70 is a complex constant. This expression, without using of the equation Q, was obt 
for the first time in P^. 

2. Let's assume in (jlj) that V = V{x, y) = 1 at —00 < x < +00, y > 0. Applying the Fi 
transformation in the variable x, we shall construct the solution ^'j (j = 1, 2, . . . , A*", . . .), ■ 
is bounded in x and goes to zero at y — > 00: 



1 

{x,y) = — / dke'^^^-^^'yA^ik). 



Here Aj{k) are some functions (functional parameters), such that Aj(k) € L2(— cxd,co 
the real case one obtains from (|19|) Aj{k) = Aj{—k) (in case of complex k we assume, 
Re\/PTT> 0). Then 

/+00 
dk e'''''-^^yAi{k)}. 
-00 

The analogue of the second relation in © yields: 

(^^^^ - ^'i^^'^):, = (^'ly^' - 
Therefore the solution ^[1] is determined by a relation Q at ^'i) = uj{^, ^i), whe: 

w(^',^'i)= / {'i'ly^ -^i^y)dx + {'i/^^i-^i^'^)dy. 

J {^0, yo) 

On the next step we obtain: 



In the particular case of Aj{k) = Cj6{k — Pj) + dj6{k — rj), with 6{.)— the Dirac delta- fun 
Cj, dj, Pj, rj £ M— constants, j = 1,2, . . . , N, . . taking (due to linearity of the equatioi 
the real part of the solution ^j we shall have: 



^j{x, y) = CjC V J cospjX + djC v ^ cosr^x. 
Thus, the dressing formula (|15|) gives: 

V[l] = l-2Aln {cie-'^P^y cospix + die-'/^y cosrix). 

Completely similarly, it is possible to obtain from (|15|) . (|24j) the explicit representation for 
though rather cumbersome. 

3. Let us consider a case, when the initial solution is determined by the Coulomb pote 
i.e. in (jlj we assume that V{x, y) = {^J + x G (— cxd, cxd), y G (— oo, oo). S( 
this equation in the polar coordinates and coming back to the cartesian ones, (with require 
that the solution is finite at x = y = 0), we obtain [TT] : 

^„(x,y) ~ e*'"-^*-^ J2,^ (2.(x2 + y^)!) , m = I, 2, . . . , N , . . . , 
where Ju{-) is the Bessel function. Then for real potential we find: 

V[l] = ^^=f - 2Aln{Re [e™-^*-^ J2,^ {2i{x^ + y^)\)] } . 

It is clear now, that already on a first step we obtain the potential which is not ame 
for separation of variables. Hence, all subsequent functions of the chain ^'[A;], V[k] obe; 
property as well. 

Potentials of the form (|2(i|) and (|28)) . and of the corresponding chains, are integrable, i.( 
equations, produced by them, have the exact solutions. 

In two following examples we consider a case gi 7^ when the equation Q allows the sin 
solutions. 

4. Let Viz = 0, i.e. is an antiholomorphic function in some domain B C C. From \ 
have: 

{\n.Vl)z = ,^ ^ at the condition, that ( ^ )z = 0. 

To solve this equation we use the formula of 9-problem (see, for example, jH]). Assi 
that 

where Co is some complex constant, we find: 



\nn{z)=Co + ^ [ [ ^ -dCAdC 



2TTiJ Jo{\n^i)^{C- z) 

Here C = C-R + ^C/i d( A dC = —2idC,B. dCi- The correctness of this expression follows froi 
well-known relation of the theory of generalized functions {d = 82) ■ 

(9(— ^ r) = S{z - zo). 

tt{z - Zo) 



5. Analogously, let = 0, i.e. is an holomorphic function in D C C. From Q we ol 



{ln^})z = ,^ ^ at the condition, that ( ^ )g = 0. 



Assuming, that 



where Ci is complex constant, from (jSHJ we have 



Inn(z) = Ci + -^T / / 7 ^ -dCAdC, 



2ni J Jo (ln^'i)^-(C-2 

d). Now we shall find the quantum- mechanical sense of the transformations (jH)). Fo 
purpose we shall consider more general, than @, equation: 

where A € C is a complex parameter. It is possible also to use here the transformation (j§ 
dressing relation (|15)1 . Let us assume that V{z,z) = zz, corresponding to the two-dimen: 
isotropic harmonic oscillator. Then the equation (|36|) can be rewritten as 

-4^^5(2, z, A) + zz'^iz, z, A) = A^'(z, z, A). 

For A = Ai = — 2 this equation has the solution ^i{z, z) = exp(^zz), and agrees witl 
V{z, z) = zz — A — Aqi. Thus we obtain the equation 

-A^,-,[l]{z, z, A) + {zz - 4(7i)^'[l](z, z, A) = (A + 4)^[l](z, z. A). 
Comparing (pTTj) and (|nH|) for (71 = 0, we find 

z, \) = ^{z, z, A + 4), 

i.e. the transformation Q, together with the requirement of covariance, acts as the quai 
mechanical creation operator of particles ^, and the quantity — 4(7i can be interpreted £ 
energy level shift (in a complex case). In this sense we have two-dimensional generalizat: 
the standard DT 

Here the equation for functional 17 follows from Q: 

^xx + ^yy ~ 22; — 2y D.y + qi^l = 0. 

This equation, which can be called the generalized Helmholtz equation, is known in the ir 
matical literature [T2I . It arises under consideration of orthogonal polynomials in two ind 
dent variables. For = 2(n + m) it has the solution 

Q{x, y) = Tn+m,m{x, v) = Hn{x)Hm{y), 

where Hi{.) are polynomials of Chebychev-Ermit, m, n = 0, 1, . . .. Here the generalized con( 
of orthogonality is 

/ / h{x, y)J^nmix, y)J^ksix, y) dx dy = 



where h{x, y) = exp{— (x^ + y^)} is the weight function, and 6ij is the Kronecker symbol 
For qi = the solution of equation (|in|) also can be factorized. Setting y) = Qi{x))i 
we have: 

- 2xO,ix - [30,1 = 0, 0,2yy " 22/1^2?/ + 13^2 = 0, 

where /3 is the "separation of variables parameter". For f3 = 2n, where n = 1, 2, . . . , the s^ 
of these equations has the solution in terms of the Chebychev-Hermit polynomials: ^2 
Hn{y)- Setting Vli{x) = V{x) = YlT=o'^>'^^^ from the first of the equations we fin 
reccurent relation: 

_ 2{k + n) 

where the numbers oq, ai remain arbitrary, and the power series for function V{x) com 
uniformly on the entire axis. 

Thus, in the dressing procedure the solution of equation ()36() is 

y) = mix, y) = Cne^'+y'v{x)Hn{y), 

where Cn is a constant. Besides that, similarly to the one-dimensional case (one-dimeni 
stationary Schrodinger equation) , the states described by the relation (|45j) are not normalr 
The detailed investigation of the creation and annihilation operators in two and more s 
dimensions was performed in the monograph 



3. A method of d - problem. 

In this Section we consider the application of the 9-problem method for the simplest c 
diffraction on half-planes (see, for example, |14j): 

Au + k^u = 0, 

and ^ 

u = -n(^)(x, y) = -e-^*^^™^^ at x < 0, y ^ ±0, 



u = 0(—;=), — iku = o(—;=) at r = \/ x"^ + y^ —>■ oo. 

^yr or yjr 

Here u = u{x, y) = n''*°*)(x, y) — u^'^\x, y), is the total potential, n^*-*) is the pot 

of a falling wave, n^*^ = exp {—ikx cos 9 — iky sin 9), A is the two-dimensional Laplace ope 
k = ki + ik2, ki, k2 > 0, < 9 < IT. 

Let us introduce the variables z and z, and let us assume u{z, z) = v{z z) exp {k-Z + 
where A;+ = —{i/2)kexp{i9), fc_ = —{i/2)kexp{—i9). Then in the z - representation the pre 
can be written in terms of function v = v{z, z) = v{z, z, k-) as 

Vzz + k+Vz + k^V2 = 0, 



v(z, z) = —1 at z + z<0, z — z^ ±iO. 



The radiation condition (|48|) gives 



f (z, z) = 0(1) at \z\ — > 00. 
It is not difficult to obtain the symmetry property from the equation (|49|1: 

v{z, z] A;_) = v{z, z; A;_, k^). 

In order to derive the integral equations we note, that for any complex-valued fui 
g{z, z) e C{B) fl Ci(D), D C C the Green formulas [S] can be used: 



9dC = - / 5c A C, f gdC= / 5C A dQ. 
Jn J Jan Jo J 

Setting 5r(C,C ^i, h) = 5i(C,C 'i, h) = w{C, C, hj l2)/{C-z), giCC h, h) = 92{CX, h, 
w{Q, C, h, l2)/iC-z), wiC, C) = wiC, C, h, h) G C(B) n C^{^). where h £ C are param 
applying to these functions the Cauchy's formula with paths OD|J — = ei, and (91D) |J |C- 
€2 accordingly, where ei, 62 > 0, and passing to limits limei^2 = 0, we obtain the in 
representations: 

w{z, z) = — — h ' ' 



27ri C ~ z Jo J 2TTi Q — z ' 
w{z, z) = - — — + ' ' 



2TTi C ~ z Jo J 27ri C ~ z 
To apply these representations, equation can be rewritten as 

dz{dz + k^)v{z, z) = -k+dz:v{z, z), d^id^ + k+)v{z, z) = -k^dzv{z, z). 
We assume, that 

Vz{z, z) = 0(1), Vz{z, z) = 0(1) at l^l — > 00. 

Choosing the path as dO = T\J Cr, where Cr is the circle of radius i?, and contour T co 
of the straight line y = itiO, at x < 0, and the half-circle of small radius at x > 0, we c 
from (jSH)-® : 



Vz{z, z) + k^v{z, z) 



dC v^ic, C) + kMC, C) f fdCAdC k+v^ic, C) 



ruCfl27ri (-z Jo J 2m C^-z 

These relations can be also simplified. Indeed, taking into account and (|57|) we have: 

/■ z;c(CC) + fc-^(C.C) ^ /■ ^c(CC)_+M(C,C) ^ Q_ 

i?-oo Jc-B 2vri C - ^ -R^oo JCn 2vrz C - ^ 

Besides that, condition gives: 



2-Ki C, — z J-oo ^ — z 

Then equations (|58j) - (|59j) yield: 



oo 



where we used notations: 



......... 1 



The equations (|6c{|) - (|64j) allow, in particular, to obtain the system of equations which cor 
a derivative of a field on a negative half-axis with derivative in the domain D. Reall 
z = X < and taking into account the boundary condition H5Ujl . equations ()63jl can be w 

as 



V- (x) - 
and equation (|HH) as 

vfix) - k+ 







6- 


X — iO 




-^ci(ei) 


6- 


X + iO 






i 6- 


- X + iO 







t^Ci A dCi ^Ci(Ci. Ci) 

27rz Ci ~ 2; 

^Ci A d(i ^Ci(Ci, Ci) 
27ri Ci ~ X 



dC2 AdC2 ^(2(^2, C2) 



27ri C2 — X 

dC2 AdC2 %(C2, C2) 



27ri £,2 — X — iO Jo J 2-iri C2 — x 

The system of equations H66() - ()69|) is closed system of singular integral equations on th 
X, to find the functions vf, v~ and vf, v^, respectively, from the values of Vz and Vz 
arbitrary point of domain B. 

These functions can be also calculated using another approach. Indeed, equation ^ 
be rewritten as 

dz + k+)vz{z, z) = -(^ dz + k^)v2{z, z). 

Then we construct the generalized function - fundamental solution of the operator ((1/4) 
idy) + k^), which satisfies the equation 

[{d^ + idy) + 4A:+]Gi(x, y, k+) = A6{x)6{y). 
Thus from the equation l|7U|l it follows: 



Similarly we can find the generalized function - fundamental solution of the operator ((1/ 
idy) + 

[{d^ - idy) + 4k.]G2{x, y,k.) = A5{x)6{y). 
From (jTUj) and (|73() we obtain: 

Vz{x,y) = - / d(,dr]G2ix-£,, y-r], k-)[d^ + id^ + 4:k+]v(;{S,, rj). 
Jo J 

The solutions Gi and G2 are connected among themselves by the symmetry relation: 

Gi{z, z, kj^) = -G2{-z, z, -k+), G2{z, z, -k-) = -Gi{-z, -z, -k-). 

The system of equations (|72j) , (|74j) is a system of integral equations of the convolution 
and it can be solved by the Fourier transformation (for generalized functions). Such calculi 
as well as the explicit expressions for functions Gi, G2 and the comparison with the resi 
jl4j will be given later. 

The approach developed above in this Section is seems to be an alternative to the tradi 
Wiener-Hopf approach. It would be interesting to generalize this method for the case of 1 
dimensional problem of diffraction. 
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